The most general fourth order theory of Gravity at low energy 

A. Stabile* 

Dipartimento di Ingegneria, Universita' del Sannio 

Corso Garibaldi, I - 80125 Benevento, Italy 

The Newtonian limit of the most general fourth order gravity is performed with metric approach 
in the Jordan frame with no gauge condition. The most general theory with fourth order differential 
equations is obtained by generalizing the f{R) term in the action with a generic function containing 
other two curvature invariants: Ricci square {RapR'' ) and Riemann square {Rap-ysR" )• The 
spherically symmetric solutions of metric tensor yet present Yukawa-like spatial behavior, but now 
one has two characteristic lengths. At Newtonian order any function of curvature invariants gives 
us the same outcome like the so-called Quadratic Lagrangian of Gravity. From Gauss - Bonnet 
invariant one have the complete interpretation of solutions and the absence of a possible third 
characteristic length linked to Riemann square contribution. From analysis of metric potentials, 
generated by point-like source, one has a constraint condition on the derivatives of / with respect 
to scalar invariants. 
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CN ■ I. INTRODUCTION 



In recent years, the effort to give a physical explanation to the today observed cosmic acceleration [i| has attracted 
O , a good amount of interest in /(i?)-gravity, where / is a generic function of Ricci scalar R, considered as a viable mech- 
^ ' anism to explain the cosmic acceleration by extending the geometric sector of field equations without the introduction 
5—1 \ of dark matter and dark energy. Other issues, od astrophysical nature, as the observed Pioneer anomaly problem [3] 
^*J' can be framed into the same approach \3\ and then, apart the cosmological dynamics, a systematic analysis of such 
theories urges at short scale and in the low energy limit. 

While it is very natural to extend the theory of General Relativity (GR) to theories with additional geometric 
K>" ■ degrees of freedom, recent attempts focused on the old idea Q of modifying the gravitational Lagrangian in a purely 
t^^ , metric framework, leading to higher-order field equations. Due to the increased complexity of the field equations in 
this framework, the main body of works dealt with some formally equivalent theories, in which a reduction of the 
order of the field equations was achieved by considering the metric and the connection as independent objects [5| . 
In addition, other authors exploited the formal relationship to scalar-tensor theories to make some statements 
C^ ' about the weak field regime fq], which was already worked out for scalar-tensor theories il\. Also a Post- Newtonian 
parameterization with metric approach in the Jordan Frame has been considered jSj] . 

In this paper, we show the theory of Gravity induced by a most general fourth order theory obtained by using 

• • all curvature invariants. Precisely we show for a generic /(AT, y, Z)-theory, where X = R, Y = RapR"^ and 

^ Z = Rap.ysR"'^'^ with Rai3 Ricci tensor and Rap-yS Riemann tensor, the modifications to standard gravitational 

k> , mechanics at Newtonian order. From the usual small velocity and weak field limit approach [9| we find the field 

\^ • equations and since the differential equations are liner we obtain a general solution of metric tensor by Green functions 

Cu I method and demonstrate that any f{X,Y, Z)-theoT'y corresponds to so-called Quadratic Lagrangian {f{X,Y) = 

aiR + a2R^ + a^RapR°^^). Initially the metric tensor is spherically symmetric and time depending, but in this limit 

the dependence is missing (we need the post-Newtonian order to fix a possible time dependence). So we recover also 

a partial outcome about the Birkhoff theorem. 

The metric potentials have two characteristic lengths depending on the value of derivatives of / with respect to 
curvature invariants and only in GR are equal. The general solutions are calculated for a point-like source and since 
the theory is linear, the gravitational potential can be obtained for any matter distribution. 

With this general approach and by adding other curvature invariants to action, this paper summarizes and gener- 
alizes the topics of previously papers |9l-[ll] . 
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II. THE NEWTONIAN LIMIT OF f{X, Y, Z)-GRAVITY 

Let us start with a general class of fourth order theories given by the action 



A= [ d-^xy/^ 



fiX,Y,Z) + XCr. 



(1) 



where / is an unspecified function of curvature invariants X, Y and Z. The term Cm is the minimally coupled 
ordinary matter contribution. In the metric approach, the field equations are obtained by varying ([1} with respect to 
g^^. We get 



H„ 



where T,, 



fxRfiv — -^5^!^ - fx:,fiiy + gfii^Ofx + 2/yi?^"i?Qi, — 2[fYR"{fi.]-„)a + 0[fY Rfii,] + [fv Ra^]'" 5^" 



'2fzRfial3jR,y'^ "' — ^[fzRtJ." ^];al3 = X T^^ 



al3 



1 ■5(v/^£„.) 



is the the energy-momentum tensor of matter, fx 



-^ fv 



rfy 



Iz 



(2) 



if 



D = :ct''^, and X — SttG^. The conventions for Ricci's tensor is i?^i/ = R'^ ^av while for the Riemann tensor is 
R'^pfjiu = r^j, „ + •■■• The affinities are the usual Christoffel's symbols of the metric: V^^ = \g^''{gaij,p+gpcr,a — gafiM)- 
The adopted signature is (H ) (see for the details [l2l)- The trace of field equations ^ is the following 



H - g-^^H^p = fxR + 2/yi?„^i?"'5 + 2fzRap-,sR"^^' - 2/ + D[3fx + fyR] + 2[(/y + 2fz)R"%o.p ^ XT (3) 

where T = T'^^ is the trace of energy-momentum tensor. 

The paradigm of Newtonian limit is starting from a develop of metric tensor (and of all additional quantities in the 
theory) with respect to dimensionless quantity v and considering only first term of tt- and ij-component of metric 
tensor g^^ (for details see [lOl)- The develop of metric tensor is as follows 



A+5*^f(t,x)+gW(i,x) + 



5,,+g(f(t,x) 



(4) 



The set of coordinates^ adopted is 



(f, x^, x^, x^). The curvature invariants X , Y, Z become 



The function / can be developed as 



Y r^ y''^'^ + yc^) - 
z - z(4) + z(6) . 



(5) 



f{X,Y,Z) ^ /(0) + /x(0)X(2) + i/x^(0)X(2)' + /x(0)xW+/y(0)y(4)+/^(0)zW + ... (6) 

and analogous relations for partial derivatives of / arc obtained. From lowest order of field equations ([2]) we have 



/(O) = 



(7) 



^ Here we use the convention c = 1. 

^ The Greek index runs between and 3; the Latin index between 1 and 3. 



Not only in /(i?)-gravity [lO; ll^l but also in /(X, F, Z)-theory a missing cosmological component in the action ([T]) 
implies that the space-time is asymptotically Minkowskian. The equations © and ([3]) at 0(2) - order become'^ 



' Hlf = /x(0)4'' - [/y(0) +4/z(0)]a42) - i^X^^^ - [/^^(O) + Mo)]AX(2) = AT^") 
i/f = /x(0)i?lf - [/y(0) + 4/z(0)]Ai?(;) + M2)x(2)5,^. + [/^^(O) + M2)]AX(2)% - /xx(0)X(2)_,^.+ 

+ [/y(0) + 4/z(0)]i?|^^„ + /y(0)i?^]-_ = 
H('' = -/x(0)X(2) _ [3/^^(0) + 2/y(0) + 2/z(0)]AX(2) = XT^^) 
where A is the Laplacian in the flat space. By introducing the quantities 



(8) 



mi 



m2 



MO) 

3/xx(0)+2/i.(0)+2/z(0) 

/x(0) 
/i'(0)+4/z(0)] 



we get three differential equations for curvature invariant X^'^' , tt- and ij-componcnt of Ricci tensor R 



(2) 



(9) 



(A - m2^)Ri 



m-2 _ mi +2m2 A 
2 tolT^ '^ 



X(2) 



m2^X rp{0) 

' fx{0) ^tt 



' {A-m2^)R, 



(2) 



mi —m,2 a2 _ 
3mi2 '^ij 



Wt2 

2 



6mi^ 



A U, 



X(2) = 



[iA-m^)X(^^ = ^-f^T(0) 



(10) 



We note that in the case of /(i?)-theory we obtained a characteristic length (mi^^) on the which the Ricci scalar 
evolves, but in f{X,Y,Z)-th.eory we have an additional characteristic length (m2^^) on the which the Ricci tensor 
evolves. The solution for curvature invariant X'^^^ in third line of (TTOl) is 



X(2)(t,x) =. 



1^ 



d3x'5i(x,x')rW(i,x') 



(11) 



where CJi(x,x') is the Green function of field operator A — mi^. The solution for ^^j , by remembering i?jj — ^Ag^^ , 
is the following 



9it\t,K) 



2n 



,,3_„52(x',x") 



d^x'd^x' 



/x(0) 6/x(0) 6 



(12) 



where Q2{x,x') is the Green function of field operator A — TO2^. The expression ([12]) is the "modified" gravitational 

potential (here we have a factor 2) for /(X, F, Z)-gravity. The solution for the gravitational potential git /'^ has a 
Yukawa-like behaviors ([10]) depending by a characteristic lengths on whose it evolves. 
The ij-component of Ricci tensor in terms of metric tensor Q is 



R 



(2) 



J2) 



.(2) 



,(2) 



.(2) 



.(2) 



tj cyyij.mm eyj irn,mj i^y jm.,mi lyjtt.ij ' 0^771771, ij 



(13) 



We used the properties: «„«•"'' = iDR and R^,°''^ ^.^f, = -Rp" ,.„ " ^Rm^ 



and if we use the harmonic gauge condition (^"^r^^ = 0) the (fT3)) becomes ([i3|) Rij \hg — ^3- 
The general solution for g^- from (llOp . in the harmonic gauge, is 



1.(2) 



1A„(2) 



2'dij,mm 2 3ij 



„(2)| 



27r 



,^3 ,,g2(x',x") 



d^x'd^x' 



.(2) 



mi" - ms" 2 _ 



3mi 



^i"j" 



( JTi'i' rrii'^ + 21712^ 
[~2 



A„,, \S^ 



Q 2 --■' r»i 



X(2)(x") (14) 



?(2) 



While if we hypothesize g,- = 2ip6ij we have i?^- = Aip Sij + {ip — 4>)^ij and the second field equation of ([TU]) 
becomes 



A^ = /d3x'g2(x,x') 



m2 
2 



mi-+2m2- ^^ \^(2)(^/) 



6?TM^ 



(15) 



^)., = ^%;^/d3x't;2(x,x')4^]-,(x') 



.(2 



Then the general solution for g\-' from ([TU]), without gauge condition and by using the first line of (|15p . is 



.(2) 



2V'(5j, 



"27r 



J |x — x'l \ 2 6mi2 y 



(16) 



and the second line of (1151) is only a constraint condition for metric potentials. In fact from its trace we have 

A(^-^) = '"^~r' /dVg2(x,x')A.,X(^)(x') 



3?7ll2 

and we can affirm that only in GR the metric potentials (j) and '0 are equals. 



(17) 



III. CONSIDERATIONS ABOUT POINT-LIKE SOLUTION 

Let us consider a point-like source with mass M. Then we have T^ (i,x) = r'^"^(i,x) = Af5(x), and if we choose 



TOi2 > and ?7i2^ > 0, the Green functions Qi become tjj(x, x') 
(fTTj) and the metric potentials (/) (fT2|) and f/' (HH) are 



47r |x— x'l 



The curvature invariant X^"^' 



(2) ^ rgmi2e-"'^N 
/x(0) |x| 



(18) 



7^ 



1 1 e~™il''l 4e~'"2|x| 



(19) 



i, 



GM 



1 1 e-»"il''l 2 e^'"^!''! 
bd " 3^d 3 Ixl 



(20) 



where rg — 2GM is the Schwarzschild radius. The modified gravitational potential by /(i?)-theory is further modified 
by the presence of functions of RapR"^ and RajS-ysR"^^^ ■ The curvature invariant X*^2) (tj-^g Rjcci scalar) presents 
a massive propagation and when f{X,Y,Z) — >■ f{R) we find the mass definition ni^ = —f'{R ~ 0)/3f"{R — 0) 
([13, Ell) and propagation mode with m2 disappear. Obviously the expressions (ITO)) and (^0]) satisfy the constraint 



* We choose a system of isotropic coordinates. 




FIG. 1: Plot of metric potential (j> (|19[). 7712 = ^?ni and mi = .1 (dotted line), mi = ^m2 and m2 = .1 (dashed line). The 
behavior of GR is shown by the solid line. The dimensionless quantity ^ runs between -^ 10 with step 2. The dimension of 
mi and m,2 is the inverse of length. We set fx{0) — 1. 




|x| 



FIG. 2: Plot of metric potential tp (|20|) . 7712 = ^m.i and m.i = .1 (dotted line), m.i = f m2 and m,2 = .1 (dashed line). The 
behavior of GR is shown by the solid line. The dimensionless quantity ^ runs between -^ 10 with step 2. The dimension of 
mi and m2 is the inverse of length. We set fx{0) ~ 1. 



condition (|17l) . In FIGs. [T] and [5] we report the spatial behavior of metric potentials for some values interval of 
parameters mi and 7712- 

The same outcome can be obtained by considering the so-called Quadratic Lagrangian C = \/—g{ai i? + 02 -R^ + 

", 1/2 
o-TiRafiR" ) where oi, 02 and 03 are constants. In this case |lll | we found two characteristic lengths 



1/2 



2(302+03) 



and the Newtonian limit of theory implied as solution the equations (IT^ and (|20p . We can affirm, then, 



the Newtonian limit of any /(X, y, Z)-theory can be reinterpreted by introducing the Quadratic Lagrangian and the 
coefficients have to satisfy the following relations 



ai - /x(0), a2^\fxx{Q)-Iz{Q) 



03 



/y(0)+ 4/^(0) 



(21) 



A first considerations about (|21l) is regarding the characteristic lengths induced by /(X, Y, Z)-theory. The second 
length m2~^ is originated from the presence, in the Lagrangian, of Ricci and Riemann tensor square, but also a 
theory containing only Ricci tensor square could show the same outcome (it is successful replacing the coefficients 
Gi of Quadratic Lagrangian or renaming the function /). Obviously the same is valid also with the Riemann tensor 
square alone. Then a such modification of theory enables a massive propagation of Ricci Tensor and, as it is well 
known in the literature, a substitution of Ricci Scalar with any function of Ricci scalar enables a massive propagation 
of Ricci scalar. We can, then, affirm that an hypothesis of Lagrangian containing any function of only Ricci scalar 



and Ricci tensor square is not restrictive and only the experimental constraints can fix the arbitrary parameters. 

A second consideration is starting from the Gauss - Bonnet invariant defined by the relation Gqb — X"^ ^ 4y + Z 
[ij]. In fact the induced field equations satisfy in four dimensions the following condition 

H'^.'' = H^u - 4^J. + H^. =. (22) 

and by substituting them at Newtonian level {Hf^ ^ —AARj-^') in the equations ^ we find the field equations (ever 
at Newtonian Level) of Quadratic Lagrangian. 

A third and last consideration is about the solutions ([T9| and (|20| . When we perform the limit in the origin |x| = 
we don't have the divergency. In fact we find 

|xKo^ 3 ' |xKo^ 3 ^ ^ 

and only if we remove in the action ([l} the dependence on the Ricci square or Riemann square we get the known 
divergence of GR. For a physical interpretation of solution (|19|) we must impose the condition mi — 4m2 < to have 
a potential well with a negative minimum in |x| = and mi < TO2 to have a negative profile of potential (see FIG. 
[1]). Then, if we suppose /x(0) > 0, we get a constraint on the derivatives of / with respect to curvature invariants 

fxx{0) + fY{0) + 2fz{0) <0 (24) 

In the case of /( /^)-g ravity (/y(0) = /z(0) = 0) we reobtain the same condition among the first and second 
derivatives of / @-[lll ■ 

IV. CONCLUSIONS 

In this paper the theory of Gravity induced by a most general fourth order theory obtained by using all curvature 
invariants has been considered. By adding these curvature invariants, this paper summarized and generalized the topics 
of previously papers [9l-[lll|. In fact for a generic f{X, Y, Z)-theory, at Newtonian level, it is successful considering only 
the so-called Quadratic Lagrangian. All contributions to field equation due by curvature invariant Riemann square 
can be expressed by other two curvature invariants (Ricci tensor square and Ricci scalar square) via Gauss-Bonnet 
invariant. 

The spherically symmetric solutions of metric tensor at C'(2)-order show a Yukawa-like dependence only by two 
characteristic lengths and not by three (because we have three curvature invariants in the action). No gauge condition 
has been considered and the solution of ij'-component of metric tensor is general. Is general also the solution of 
tt-component of metric tensor, since it is gauge free (but only at Newtonian order). 

Furthermore, generally it has been shown that for a f{X, Y, Z)-gravity, but the same is valid also for /(i?)-gravity, 
the metric potentials are not equal. Only in the limit f ^ R we obtain the outcome of GR. This aspect with the 
consequences of Birkhoff and Gauss theorem are the principal differences between a fourth order gravity and GR. 

The general solutions are calculated for a point-like source and since the theory is linear, the gravitational potential 
can be obtained for any matter distribution. The metric potentials don't have the divergency in |x| = for a point- 
like source, and by requiring a right physical interpretation of solution we get a constraint on the derivatives of / with 
respect to curvature invariants. The constraint condition is compatible with respect to one obtained for /(i?)-gravity. 
Besides a such class of theory have free parameters and only the experimental evidence can fix them. 

V. ACKNOWLEDGES 

The author would like to thank prof. Salvatore Capozziello for his useful discussions and suggestions. 



[1] A.G. Riess et al. Astron. J. 116, 1009 (1998) 

S. Perlmutter et al. Astrophys. J. 517, 565 (1999) 



D.N. Spergel et al. Astrophys. J. Suppl. 148, 175 (2003) 

A.G. Riess et al. Astrophys. J. 607, 665 (2004) 

S. Cole et al, Mon. Not. Roy. Astron. See. 362, 505 (2005) 

S. Perlmutter et al. Astron. Astrophys. 447, 31 (2006) 

D.N. Spergel et al. Astrophys. J. Suppl. 170, 377 (2007) 
[2] J.D. Anderson et al. Phys. Rev. Lett. 81, 2858 (1998) 

J.D. Anderson et al. Phys. Rev. D 65, 082004 (2002) 
[3] O. Bertolami, C.G. Bohmer, T. Harko, F.S.N. Lobo, Phys. Rev. D 75, 104016 (2007) 
[4] H. Weyl, Raum-Zeit-Materie, Springer Berhn (1921) 

A.S. Eddington, The mathematical theory of relativity, Cambridge University Press London (1924) 

C. Lanczos, Z. Phys. 73, 147 (1931) 

H.A. Buchdahl, Nuovo Cimento 23, 141 (1962) 

G.V. Bicknell, Journ. phys. A 7, 1061 (1974) 

P. Havas, Gen. Rel. Grav. 8, 631 (1977) 

K.S. Stelle, Gen. Rel. Grav. 9, 353 (1977) 
[5] G. Magnano, M. Ferraris, M. Francavigha, Gen. Rel. Grav. 19, 465 (1987) 

G. AUemandi, A. Borowiec, M. Francavigha, Phys. Rev. D 70, 103503 (2004) 

M. Amarzguioui, O. Elgaroy, D.F. Mota, T. Muhamaki, Astron. Astrophys. 454, 707 (2006) 
[6] G. J. Olmo, Phys. Rev. D 72, 083505 (2005) 

G. J. Olmo, Phys. Rev. Lett. 95, 261102 (2005) 

G. J. Olmo, Phys. Rev. D 75, 023511 (2007) 
[7] T. Damour, G. Esposito-Farese, Class. Quant. Grav. 9, 2093 (1992) 
[8] T. Clifton, Phys. Rev. D 77, 024041 (2008) 

[9] S. Capozziello, A. Stabile, A. Troisi, Phys. Rev. D 76, 104019 (2007) 
[10] S. Capozziello, A. Stabile, A. Troisi, Modern physics letters A 24, 659 (2009) 

A. Stabile, arXiv: gr-qc/1004.1973v2 
[11] S. Capozziello, A. Stabile, Class. Quant. Grav. 26, 085019 (2009) 
[12] Landau Lev D., Lifsits E. M, Theoretical physics vol. II 
[13] S. Capozziello, A. Stabile, A. Troisi, Class. Quant. Grav. 25, 085004 (2008) 
[14] B.S. de Witt, Dynamical theory of groups and fields, Gordon and Breach, New York (1965) 



